Limits on Spacetime Foam 



Wayne A. Christiansen and Y. Jack Ng 
Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599 

David J. E. Floyd 
AAO/OCIW Magellan Fellow. Current address: School of Physics, 
University of Melbourne, Victoria 3010, Australia 

Eric S. Perlman 

Physics and Space Sciences Department, Florida Institute of Technology, Melbourne, FL 32901 

Plausibly spacetime is "foamy" on small distance scales, due to quantum fluctuations. We elab- 
orate on the proposal to detect spacetime foam by looking for seeing disks in the images of distant 
quasars and AGNs. This is a null test in the sense that the continued presence of unresolved "point" 
sources at the milli-arc second level in samples of distant compact sources puts severe constraints 
on theories of quantized spacetime foam at the Planckian level. We discuss the geometry of foamy 
spacetime, and the appropriate distance measure for calculating the expected angular broadening. 
We then deal with recent data and the constraints they put on spacetime foam models. While time 
lags from distant pulsed sources such as GRBs have been posited as a possible test of spacetime 
foam models, we demonstrate that the time-lag effect is rather smaller than has been calculated, due 
to the equal probability of positive and negative fluctuations in the speed of light inherent in such 
models. Thus far, images of high-redshift quasars from the Hubble Ultra-Deep Field (UDF) provide 
the most stringent test of spacetime foam theories. While random walk models (a = 1/2) have 
already been ruled out, the holographic (a = 2/3) model remains viable. Here a ~ 1 parametrizes 
the different spacetime foam models according to which the fluctuation of a distance I is given by 
~ l l ~ a l P with lp being the Planck length. Indeed, we see a slight wavelength-dependent blurring in 
the UDF images selected for this study. Using existing data in the Hubble Space Telescope (HST) 
archive we find it is impossible to rule out the a = 2/3 model, but exclude all models with a < 0.65. 
By comparison, current GRB time lag observations only exclude models with a < 0.3. 



I. INTRODUCTION 

Even at the minute scales of distance and duration 
examined with increasingly discriminating instruments, 
spacetime still appears to be smooth and structureless. 
But, like everything else, (plausibly) spacetime is sub- 
ject to quantum fluctuations. So we expect that space- 
time, probed at a small enough scale, will appear com- 
plicated — something akin in complexity to a turbulent 
froth that [45[ has dubbed "quantum foam," also known 
as "spacetime foam." But how large are the fluctua- 
tions in the fabric of spacetime? How small a scale is 
small enough so that the foaminess of spacetime mani- 
fests itself? To quantify the problem, let us recall that, if 
spacetime indeed undergoes quantum fluctuations, there 
will be an intrinsic limitation to the accuracy with which 
one can measure a macroscopic distance, for that dis- 
tance fluctuates. Denoting the fluctuation of a distance I 
by 51, on general grounds, we expect 61 > Nl 1 " ^, (23j 
where N is a numerical factor ~ 1 and lp = y KG / c? 
is the Planck length, the characteristic length scale in 
quantum gravity, and we have denoted the Planck con- 
stant, gravitational constant and the speed of light by 
h, G and c respectively. It is important to note that 51 
cannot be defined without reference to the macroscopic 
distance, I (i.e., 81 is not defined locally). The parameter 
a < 1 specifics the different spacetime foam models. We 
note that smaller values of a necessarily lead to larger 



distance fluctuations and, hence, such models are easier 
to test observationally. 

Because the Planck length (~ 10 -33 cm) is so incredi- 
bly small, we need an astronomically large distance I for 
its fluctuation 51 to be detectable. Even so, measurement 
of 51 is not trivial. It must be calibrated against a known 
length standard. The most obvious calibrator is the 
wavelength A, of the light received from a distant source. 
Thus, in principle, distance fluctuations ±51, imply phase 
fluctuations ±A</> = ±2tt51/X (c.f. [H, Elf). More 
recently, Christiansen et al. Q, henceforth CNvD, sug- 
gested that fluctuations in the direction of the local wave- 
vector, ±5ij) = ±A0/(2tt) = ±51/ X could possibly be de- 
tected as "halos" in images of distant "point sources". 
However, it must be noted that using the wavelength of 
received light to bootstrap a calibration for measuring 
51, necessarily involves interferometry and/or wave op- 
tics such as Strehl analysis (see section ITV|) . The point 
is that, due to quantum foam- induced fluctuations in the 
phase velocity of an incoming light wave from a distant 
point source, the wave front itself develops a small scale 
"cloud of uncertainty" equivalent to a "foamy" structure. 
This results in the wave vector, upon detection, acquiring 
a jitter in direction with an angular spread of the order 
of 5yj. In effect, spacetime foam creates a "seeing disk" 
whose angular diameter is 
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For a telescope or interferometer with baseline length D, 
this means that dispersion (on the order of Sip in the 
normal to the wave front) will be recorded as a spread 
in the angular size of a distant point source, causing a 
reduction in the Strehl ratio, and/or the fringe visibil- 
ity when dip ~ X/D for a diffraction limited telescope. 
Thus, in principle, for arbitrarily large distances space- 
time foam sets a lower limit on the observable angular 
size of a source at a given wavelength A. Furthermore, the 
disappearance of "point sources" will be strongly wave- 
length dependent happening first at short wavelengths. 
In fact, as we show below, the frequency dependence of 
blurring currently offers the best chance to detect the 
effects of spacetime foam. 

As discussed by CNvD, there are three major theoreti- 
cal models for specifying the cumulative effects of space- 
time foam: 

• The random-walk model in which the spacetime 
foam effects grow like a one-dimensional random 
walk, corresponding to a = 1/2, N = 1. 

• The holographic model [IJ Hfl H3 which is con- 
sistent with the holographic principle [U [42[ 
where the information content in any three- 
dimensional region of space can be encoded on a 
two-dimensional surface surrounding the region of 
interest, corresponding to a — 2/3 and 0.97 < N < 
1.9 as discussed in section HT1 

• An anti-correlation model [20j in which there are 
no cumulative effects, so the distance fluctuation 
remains simply the Planck length, corresponding 
to a= 1,N= 1. 

The random walk model can be ruled out on the ba- 
sis of HST observations of relatively nearby quasars (dis- 
tance ~ 1 Gpc) according to CNvD. The predicted space- 
time foam induced seeing disks (~ few xlO -4 radians) 
clearly exceed the HST resolution even in the IR wave- 
band at 800 nm (3.3 x 10~ 7 radians). Although the holo- 
graphic model's seeing disk (~ 10 -8 radians) is still an 
order of magnitude smaller than the nominal IR HST res- 
olution, at larger distances (redshifts z > 1) and shorter 
wavelengths HST may have sufficient resolution to test 
this model also. On the other hand the anti-correlation 
model with Sip ~ Ip/X will likely remain untestable by 
astronomical means since the effect is vanishingly small 
and independent of distance. 

Testing models of spacetime foam requires large dis- 
tances to maximize the effect which, therefore, requires 
that cosmological effects must be taken into account. 
The fundamental cosmological distance measure is the 
co-moving distance (c.f. 10]). However, other distance 
measures may be used, depending on circumstances. For 
example, the luminosity distance is commonly used to 
interpret observations of high redshift objects (as used, 
for example, by [39| . in his study of spacetime foam an- 
gular broadening). However luminosity distance is inap- 
propriate for studying spacetime foam induced angular 



broadening. The reason is that the luminosity distance 
has a built-in energy dependence because it is based on 
the flux detected by a distant observer. Therefore, the 
redshift dependence of the luminosity distance explicitly 
includes a correction for the diminished energy received 
from distant sources. That enhances the expected an- 
gular broadening (c.f. Eq. (HJ). However, the angular 
broadening caused by spacetime foam does not explicitly 
depend on flux. Another distance measure, the angular 
diameter distance, may be appropriate for observations 
involving lensing or scattering. However, spacetime foam 
angular uncertainties are not associated with scattering 
or lensing (see section [IV] for detailed discussion). Quan- 
tum foam is, by definition, a property of the fabric of 
spacetime itself, and not an extrinsic property measured 
against the background geometry. Therefore the appro- 
priate cosmological distance measure is the fundamental 
co-moving distance which we shall use in our subsequent 
calculations. 

In addition, other factors also make the search for ob- 
servable effects more complicated for the following rea- 
sons: 

1. Sensitivity: More distant sources are fainter. This 
means that, as far as interferomctry is concerned, 
the lack of observed fringes may simply be due to 
the lack of sufficient flux rather than the possibility 
that the instrument has resolved a spacetime foam 
generated halo. This is probably the most serious 
limitation on experiments to search for spacetime 
foam. As discussed above, the effect is cumulative, 
requiring large distances to appear. But large dis- 
tances inevitably mean faint sources. 

2. Atmospheric turbulence: Integration times must be 
shorter than atmospheric fluctuation times. Thus 
again what one needs is to find a bright but suf- 
ficiently distant external source, for then the inte- 
gration times may be short enough to allow one to 
remove or at least minimize the atmospheric effects. 
Of course, the latter problem does not intrude for 
space-based observatories such as HST or Chandra, 
nor for an interferometer. 

3. Masking: This may originate under two distinct 
circumstances: 

(a) The source may have an intrinsic structure 
comparable in angular size to the expected 
spacetime foam halo; 

(b) Any competing scattering mechanism in the 
intervening medium between source and ob- 
server (e.g., gravitational wave scattering of 
light, micro-lensing by machos or dark matter, 
etc.) 

Concerning (a) , distant sources may have extended struc- 
tures. Thus if fringes are not seen, it may be an indica- 
tion that the source quasar has an intrinsic core structure 
whose angular size is comparable to, or larger than, the 
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spacetime foam halo. Since we do not know the physi- 
cal scales of all components of AGN and quasars we may 
find that spacetime foam effects are masked by intrinsic 
structures. However, since the angular size of intrinsic 
structures decreases like / _1 (in a flat universe) whereas 
spacetime foam halos increase as Z 1_Q , statistical samples 
can, in principle, separate the two effects. 

The bottom line is that a number of modern tele- 
scopes including the Very Large Telescope VLTI, HST 
and Chandra are on the verge of testing theories of space- 
time foam (see, e.g., [3{|), providing a sufficiently large 
sample of sufficiently bright compact sources can be de- 
veloped. This is where future large aperture telescopes 
will prove to be crucial. The test is simply a question of 
the detection of a reduction in Strehl ratio or correspond- 
ingly a reduction in interferometer fringe amplitude. It 
is not a question of mapping the structure of the pre- 
dicted seeing disk. Knowing the real sensitivity limits of 
the telescopes involved will enable us to search the cat- 
alogs to see if we can find some quasars that exceed the 
sensitivity limits. 

One of the purposes of this paper is to give a more 
quantitative analysis of the proposal just described. For 
that, the order of magnitude estimate for SI available so 
far is not enough. In the next section, we will calculate 
the numerical factor in the SI expression. 

A recent observation [l| of a noticeable spread in ar- 
rival times for high energy gamma rays from distant 
gamma ray bursts has prompted the suggestion that it is 
unlikely that fuzziness of radio or optical images of dis- 
tant extragalactic sources would be observed. The latter 
observation is exactly what this paper is about, thus we 
need to show that we disagree with the suggestion. Our 
detailed argument (based on previous work [251 ]) is given 
in section ITTT1 

In section UVl we present a more detailed discussion of 
wave vector uncertainty in spacetime foam models. We 
discuss specifics of the necessary null tests for spacetime 
foam models as well as the reasons why spacetime foam 
halos cannot be interpreted as resulting from multiple 
scattering events. 

So far, in our attempts to detect spacetime foam, 
only the HST images of a few quasars/ AGNs [3(3, HH 
have been analyzed [(| [ID, [24| [3{J. In section [V] we 
show the size of expected spacetime foam broadening 
for various telescopes in a variety of wavebands. We 
then concentrate on the optical regime. Noting that the 
VLTI will not be in full operation for some time, we 
concentrate our attention in this paper on an analysis 
of quasar images in the HST archive. This is done 
in section [VlJ We discuss results and conclusions in 
section I VIII 



II. DISTANCE FLUCTUATION EXPRESSION 

In this section, we will calculate the numerical factor N 
in the distance fluctuation expression for 81 > Nlpl 1 ^ 01 . 

(a) For the random walk model with step length lp, tak- 
ing M steps with a probability p (q) of taking a 
forward (backward) step, the mean square displace- 
ment is given by (SI) 2 = AMpqlp. For the distance 
I, the number of steps is M — l/lp. Therefore, with 
an unbiased random walk, p = q = 1/2, we get an 

root-mean-square displacement of SI = l 1 ^ 2 ^ 2 , so 
N=l. 

(b) Obviously for the anti-correlation model, since SI is 
always equal to lp, again N = 1. 

(c) For the holographic model, we now have an explicit 
calculation which sets narrow limits for N. We use 
four different methods and their consistency with 
each other (to within a factor of two or so, as we 
will show) bodes well for the robustness of our re- 
sult. The ingredients behind all the approaches are 
black hole physics and quantum mechanics. 

A. Quantum computation 

This method 0, [lj| relies on the fact that quantum 
fluctuations of spacetime manifest themselves in the 
form of uncertainties in the geometry of spacetime. 
Hence the structure of spacetime foam can be inferred 
from the accuracy with which we can measure that 
geometry. Let us consider mapping out the geometry 
of spacetime for a spherical volume of radius I over the 
amount of time T = 21 /c it takes light to cross the vol- 
ume. One way to do this is to fill the space with clocks, 
exchanging signals with other clocks and measuring the 
signals' times of arrival. This process of mapping the 
geometry of spacetime is a kind of computation, in which 
distances are gauged by transmitting and processing 
information. The total number of operations, including 
the ticks of the clocks and the measurements of signals, 
is bounded by the Margolus-Levitin theorem [18| in 
quantum computation, which stipulates that the rate of 
operations for any computer cannot exceed the amount 
of energy E that is available for computation divided 
by irh/2. A total mass M of clocks then yields, via 
the Margolus-Levitin theorem, the bound on the total 
number of operations given by (2Mc 2 /irh) x 21/ c. But 
to prevent black hole formation, M must be less than 
lc 2 /2G. Together, these two limits imply that the total 
number of operations that can occur in a spatial volume 
of radius / for a time period 21 /c is no greater than 
2(l/lp) 2 /tt. To maximize spatial resolution, each clock 
must tick only once during the entire time period. If we 
regard the operations partitioning the spacetime volume 
into "cells" , then on the average each cell occupies a spa- 
tial volume no less than (4tt1 3 /3)/(21 2 /t:1 2 p ) = 2tt 2 U 2 p /3, 
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yielding an average separation between neighboring cells 
no less than (2ir 2 /S) 1 / 3 / 1 / 3 /^/ 3 . This spatial separation 
is interpreted as the average minimum uncertainty in 
the measurement of a distance I, that is, 81 > 1.92 1//3 Zp ; 
thus, N = 1.9. 



B. Heisenberg's Uncertainty Principle 

This method involves an estimate of the maximum 
number of (massless) particles that can be put inside 
a spherical region of radius I. Heisenberg's uncertainty 
principle can be invoked to give the minimum mo- 
mentum of each particle to be | x (7j/2Z), or in other 
words, the minimum energy of each particle is no less 
than i x hc/l. To prevent black hole formation, the 
total energy is bounded by i x lc 4 /G. Thus the total 
number of particles is no greater than 2(l/lp) 2 . This 
yields SI > (?f) 1 / 3 x l 1 / 3 ^ 3 = 1.3l^ 3 l P /3 ; thus, N = 1.3. 



a beam of light. For the uncertainty in distance not to 
exceed 81, the clock must tick off time fast enough so 
that d/c < 81 /c. But d must be larger than twice the 
Schwarzschild radius 2Gm/c 2 . These two requirements 
imply 81 > AGm/c 2 which combines with the requirement 
from quantum mechanics to yield (SI) 3 > All 2 p (indepen- 
dent of the mass m of the clock), or SI > 1.6Z 1//3 ?p , 
thus, N = 1.6. A remark is in order: taking into ac- 
count also the jiggling of the mirror's position can only 
enhance SI. As another remark, we note that the prob- 
lem of distance measurement has been discussed in [7| in 
non-critical string theory. 

Thus, all four different methods yield very similar 
results. (Part of the relatively small discrepancy can 
be traced to the mismatch between the geometries of 
spheres and cubes.) This is not surprising as our result 
depends not on any particular theory of quantum gravity, 
but only on the quantum nature of measurements. In the 
discussion below, we will use the preceding arguments to 
set the allowable range for N in the holographic model 
as : 0.97 < N < 1.9. 



C. Holographic Principle 

This method [H, HH is based on the upper bound 
on the number / of degrees of freedom that can be 
put in a spherical volume of radius I. According to 
the holographic principle, [U |42[ the entropy S in the 
sphere is bounded by S/ks < \A-k1 2 /l 2 P , where fc# is 
the Boltzman constant. Recalling that e s l kB = 2 1 , 
we get / < ir(l/lp) 2 /ln2. The average separation 
between neighboring degrees of freedom then yields 
Si > (4te2)V3ji/3 Z 2/3 = o.g^ 1 / 3 ^ 73 ; thus, N = 0.97. 



D. The Wigner-Salecker Gedanken Experiment 

In this method [H, H3, [M H?} we conduct an ideal 
thought experiment to measure a distance I between a 
clock and a mirror, By sending a light signal from the 
clock to the mirror in a timing experiment, we can de- 
termine I. The clock's and the mirror's positions jiggle 
according to Heisenberg's uncertainty principle, resulting 
in an uncertainty SI in the measurement of I. 

Let us concentrate on the fluctuation of the distance 
I due to the jiggling of the clock's position. If the clock 
of mass m has a linear spread SI (and hence an uncer- 
tainty of speed, via Heisenberg's uncertainty principe, 
given by h/(2mSl)) when the light signal leaves the clock, 
then its position spread grows to 81 + K 2 x g , — when the 
light signal returns to the clock, with the minimum at 
51 = hi/ mc. Thus quantum mechanics alone would sug- 
gest using a massive clock. Now consider the clock to 
be light-clock consisting of a spherical cavity of diameter 
d, surrounded by a mirror wall between which bounces 



III. SPREAD IN ARRIVAL TIMES FROM 
DISTANT PULSED SOURCES 

In this section we will discuss the technique of using the 
spread in arrival times of photons as a possible technique 
for detecting spacetime foam. We are motivated by the 
interesting detection of a minimal spread in the arrival 
times of high energy photons from distant GRB reported 
by [l| ■ While useful in putting a limit on the variation of 
the speed of light of a definite sign, this technique is far 
less useful than the measured angular size in constraining 
the degree of fuzziness of spacetime in the spacetime foam 
models that we consider in this paper. The reason is that 
spacetime foam models predict that the speed of light 
fluctuates (i.e., the fluctuation takes on ± sign with equal 
probability): at one instant a particular photon is faster 
than the average of the other photons, but at the next 
instant it is slower. We will show that the cumulative 
effect of the fluctuations is too small to be detectable , 
even for high-energy gamma rays that have travelled ~ 
7 billion light-years halfway across the Universe [l|, for 
most of the spacetime foam models. 

First we need to examine the cumulative effects [24[ 
of spacetime fluctuations over a large distance. Consider 
a distance I, and divide it into l/X equal parts each of 
which has length A (A can be as small as lp). If we start 
with SX from each part, the question is how do the l/X 
parts add up to SI for the whole distance I. In other 
words, we want to find the cumulative factor C defined 
by SI = CSX, Since 61 ~ l^lp and SX ~ X^lp, the 

result isC= (j) 1 Note that the cumulative factor 
is not linear in (l/X), i.e., |^ ^ j. (In fact, it is much 
smaller than l/X). The reason for this is obvious: the 
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<5A's from the l/X parts in I do not add coherently. 

Next we have to find how much the speed of light fluc- 
tuates. We proceed as follows: Just as there are uncer- 
tainties in spacetime measurements, there are also un- 
certainties in energy-momentum measurements due to 
spacetime foam effects. Thus there is a limit to how 
accurately we can measure and know the energy and mo- 
mentum of a system (26l . |27| . Imagine sending a particle 
of momentum p to probe a certain structure of spatial 
extent I so that p ~ K/l. It follows that Sp ~ (h/l 2 )8l. 
Spacetime fluctuations 81 > l(lp/l) a can now be used to 
give 

Sp~ p (jL.Y. ( 2 ) 
\mpcj 

The corresponding statement for energy uncertainties is 

(3) 

Here mp and -Ep = mpc 2 denote the Planck mass and 
energy respectively. These energy-momentum uncertain- 
ties modify the dispersion relation for the photons to 
read: 

EP^w+dp^y, (4) 

where e ~ ±1. Thus the speed of light v — ^ fluctuates 
by 

8v ~ 2ec{E/E P ) a . (5) 

We emphasize that all the fluctuations take on ± sign 
with equal probability (like a Gaussian distribution about 
c). 

It follows from the above discussion that for pho- 
tons emitted simultaneously from a distant source com- 
ing towards our detector, we would expect an energy- 
dependent spread in their arrival times. To maximize 
the spread in arrival times, we should look for energetic 
photons from distant sources. So the idea is to look for 
a noticeable spread in arrival times for such high energy 
gamma rays from distant gamma ray bursts. This pro- 
posal was first made by Q in another context. 

To underscore the importance of using the correct cu- 
mulative factor to estimate the spacetime foam effect, 
let us first proceed in a naive manner. At first sight, 
the fluctuating speed of light would seem to yield an 
energy-dependent spread in the arrival times of pho- 
tons of energy E given by St ~ 5v{l/c 2 ) ~ t(E/Ep) a , 
where t = l/c is the average overall time of travel from 
the photon source (distance I away). But these results 
for the spread of arrival times of photons are not cor- 
rect, because we have inadvertently used l/X ~ Et/H 
as the cumulative factor instead of the correct factor 
(l/X) 1 ^" ~ {Et/K) l ~ a . Using the correct cumulative fac- 
tor, we get a much smaller 

6t~t 1 - a tf.~8l/c (6) 



for the spread in arrival time of the photons [25], inde- 
pendent of energy E (or photon wavelength, A). Here 
tp ~ I0 -44 sec is the minuscule Planck time. Thus the 
result is that the time-of-flight differences increase only 
with the (1— a) -power of the average overall time of travel 
t = l/c from the gamma ray bursts to our detector, lead- 
ing to a time spread too small to be detectable (except 
for the uninteresting range of a close to 0). 

The new Fermi Gamma-ray Space Telescope results [l[ 
of St < 1 s for t ~ 7 billion years, rule out only spacetime 
foam models with a < 0.3. We note that the result St > 
t 1-a tp is to be expected if we recall that the spread of 
arrival times can be traced to uncertainty in the distance 
that the photons have travelled from the distant source 
to our telescopes. In other words, gamma ray arrival time 
dispersions ~ 1 s from sources at distances ~ 7x I 9 ly 
only reject spacetime foam models in which the effects 
almost add coherently (0 < a < 0.3) but do not test 
the major models discussed above (i.e., a = 1/2, 2/3, 1 ). 
For example the holographic model predicts an energy 
independent dispersion of arrival times ~ 2.5 x 1 0~ 24 s. 

Finally, in the following section on spacetime foam in- 
duced uncertainties in wave vector direction we show an- 
gular uncertainties in the wave vector direction are not 
the result of random multiple scattering events which 
could, in principle, lead to much larger uncertainty in 
arrival times. 



IV. ANGULAR FLUCTUATIONS: 
MEASUREMENT 

At the outset it must be re-emphasized that the test 
for spacetime foam effects is a null test. A theoretical 
model for spacetime foam is disproved if images of a dis- 
tant point source do not exhibit the blurring predicted 
by theoretical spacetime foam models. 

As discussed in the Introduction, the fundamental un- 
certainties caused by spacetime foam are spatial, not an- 
gular. Strictly speaking, the models specify the uncer- 
tainty ±81, in distance between a source and observer 
along the line of sight. This is because SI is defined by 
the uncertainty in the distance measured by light travel 
times. Straightforward extensions of the models hypoth- 
esize that SI is isotropic in the sense that the uncertainty 
perpendicular to the line of sight is equal to the line of 
sight uncertainty (cf. CNvD, 2006). Of course, there 
is also a corresponding uncertainty in the transit time 
for light from source to observer, St ~ 51 /c, which was 
discussed in section III. Furthermore, since the globally 
averaged wavefront is effectively spherical it implies that 
globally averaged photon trajectories deviate from the 
direct line of sight by an angle less than or equal to 51/1 
. A direct consequence of this limitation on the globally 
averaged angular deflection of photon trajectories is that 
the theoretically expected blurring of distant images due 
to spacetime foam is not the result of a random walk 
of small angle photon scatterings along the line of sight, 
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since the uncertainties in the derived directions of the 
local wave vectors must result in the same spatial un- 
certainty, 51 (no matter how many wave crests pass the 
observer's location). For example, in the "thin screen ap- 
proximation" , the accumulated transverse path of mul- 
tiply scattered photons is approximated as (Stp)l » 81. 
This would lead to expected time lags, 8ip(l/c) » Sl/c, 
in conflict with the basic premises for space time foam 
models. 

The reasons why a straight- forward application of stan- 
dard small angle scattering theory (e.g., gravitational 
micro-lensing) is not applicable are two-fold: 
(1) the "scattering centers" for spacetime foam would 
be virtual Planck masses that exist for only a Planck 
time (~ 10~ 44 sec) which is much less than the period 
(~ 1CT 28 sec) of even the highest energy (TeV) detectable 
gamma rays. 

(2) the "density" of such virtual scattering centers can 
only be inferred from nonlocal theories such as the holo- 
graphic principle. . 

However, according to theory, although the wavefront 
remains approximately spherical its precise location and 
structure is best characterized as a " cloud" of uncertainty 
whose scale is 61. These uncertainties are not realized as 
fluctuations until physical measurements are performed. 
This is consistent with the quantum physics of these mod- 
els with their reliance on the quantum Uncertainty Prin- 
ciple. 

As the light wave encounters a "detector" (e.g. a tele- 
scope) the spatial uncertainty shows up as a random 
series of fluctuations of scale < (±)5x along and per- 
pendicular to the line of sight causing fluctuations (i.e. 
"jitter") in the apparent angular direction of the wave 
vector recorded by the detector. However, the cloud of 
uncertainty associated with these spacetime foam models 
requires that the spatial uncertainties in the wave, from 
crest to crest, do not continuously connect. This cloud 
of uncertainties is chaotic and fluctuations are not cor- 
related from crest to crest. That is, you cannot trace 
caustics from one crest to another. 

Furthermore, it must be understood that a necessary 
condition for measurement of angular blurring requires 
that 61 << A, which corresponds to the limit for quan- 
tum coherence. As the electromagnetic wave from the 
extragalactic source propagates, the phase uncertainty 
A<p ~ 2tt81/X gets larger. Once A<p ~ 2tt (corresponding 
to 81 <~ A) quantum coherence is lost. At this point, while 
the amplitude of the electromagnetic wave remains in- 
tact, all phase information is lost. The distance at which 
spacetime foam induced de-coherence sets in depends on 
the model. For the "a-models," coherence is lost sooner 
when a is smaller. 

In the realm where imaging is possible (81 << A), de- 
tecting spacetime foam effects via "blurring" must in- 
volve a form of triangulation by comparing 81 with some 
well-defined localized length scale, B, that would allow 
for the measurement of an angular uncertainty defined 
by the dimensionless ratio, Sip ~ 81/ B. Possible length 



scales defining B that may be relevant are: 1) the di- 
ameter of the telescope, D, or 2) the wavelength of the 
radiation, A. In both cases the angular uncertainty would 
represent the theoretically predicted apparent broaden- 
ing of the image of a point source of radiation caused by 
spacetime foam as recorded by a telescope plus detector 
GEOMETRIC OPTICS: 

If the telescope/detector does not record phase infor- 
mation, then angular broadening at the focal plane will 
be determined by ray-tracing. Therefore, the spatial un- 
certainty, 81, in the location of a ray on the primary 
(whether it is a lens or a mirror), will result in an an- 
gular uncertainty, at the focal plane of, 

8lp\Geometric Optics ~ 8lp\ min ~ 281/ D « X/D. (7) 

(Of course the focal ratio of the telescope also plays 
a role in the above equation, but the key parameter is 
the size of the primary.) When combined with the re- 
quirement for quantum coherence discussed above (i.e., 
81 << A), eq. (7) virtually guarantees that spacetime 
foam effects will be undetectable by techniques that are 
phase insensitive (i.e. geometric optics). 

WAVE OPTICS: 

Clearly the largest measurable angular effect would be 
associated with using the shortest length standard, i.e. 
the wavelength. As stated previously, the only way to 
accomplish such a measurement is via interference tech- 
niques such as diffraction limited imaging and/or inter- 
ferometry. In this case, since 81 = Nl^- a H% the angular 
uncertainty (i.e. the blurring angle) is, 

S^\wave Optics = <*VW = Sl/X = N(l/X)^ (l P /X) a . 

(8) 

This relation is identical in form to the relation defining 
81, with dimensionless ratios, (I, lp)/X, replacing the pre- 
viously defined distances. Note the so-called "accumula- 
tion factor, (l/X)^ a \ plays the same role as discussed 
previously in the growth of the distance uncertainty and 
the time uncertainty discussed in section III. The mag- 
nitude of Sip as given in Eq. (8) is consistent with our 
assumption of isotropic fluctuations which implies com- 
parable sizes of the wave-vector fluctuations perpendic- 
ular (Skp er p) to and along (8k par ) the line of sight (cf. 
CNvD, 2006). To wit, the size of the angular spread Sip 
is given by Yi5k perp / k where £ stands for summation over 
the l/X intervals each of length A. Since Sk perp ~ 8k par 
according to the assumption of isotropic fluctuations and 
kpar ~ k = 2n/X, we have Sip ~ T,Sk/k = £<5A/A ~ Sl/X. 

Thus, interference techniques offer the best (perhaps 
only) method for testing for the presence of spacetime 
foam because the observed wavelength offers the short- 
est line of sight baseline calibrator for measuring the lo- 
cal angular uncertainty in wave vector direction and in- 
terference is the only way to exploit the wavelength as 
a calibrator. That is, the expected blurring of diffrac- 
tion limited images is a consequence of interference cf- 
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fects resulting from wavefront "distortions" (or "mea- 
surement errors") corresponding to ratios of positional 
uncertainties of order ±61 along, and perpendicular, to 
the line of sight and the corresponding angular uncer- 
tainties ±Sl /A. In particular, the angular uncertainty 
corresponds to the "wavefront errors" that are the ba- 
sis for calculating Strchl ratios in wave-optics. Since the 
Strehl ratio is designed to provide a quantitative measure 
of the departures from diffraction limited optics, the tech- 
nique is well suited for providing the proposed null test 
of spacetime foam models. Simply put, if distant point 
sources do not show reductions in Strehl ratios predicted 
by theoretical models, the model must be rejected. The 
application of this technique is discussed in detail in sec- 
tion VI below. 

We end this section with a summary of what we ex- 
pect a telescope would observe at different locations. As 
the electromagnetic wave from the extragalactic source 
propagates, the phase uncertainty Acj> ~ 2nSl/X gets 
larger. So long as A<j> << X/D (in other words, so 
long as SI << X 2 /D), normal interference pattern is ob- 
served. Farther from the source, as A<f> ~ X/D (i.e., 
61 ~ A 2 /D), we expect the interference pattern to be re- 
placed by a halo structure. But once A</> reaches ~ 2ir 
(corresponding to 81 ~ A) , quantum coherence is lost, 
and the images will become all blurred. At this point, 
while the amplitude of the electromagnetic wave remains 
intact, all phase information is lost. How far from the 
source these different stages (from interference pattern 
to halo structure to all blurred images) are encountered, 
of course, depends on the different spacetime foam mod- 
els. For example, for the AGN PKS 1413+135 [H, the 
random-walk model predicts that Acf> ~ 10 x 2tt, while the 
holographic model yields Atfi ~ 10~ 9 x 2ir for A « 1.6/1 m 
@. Thus, according to the random- walk model, A(f> has 
reached the stage beyond the quantum coherence limit 
and hence the image should be all blurred. However, even 
supplemented with the assumption of isotropic fluctua- 
tions (viz., 6k perp ~ 8k par ) which henceforth we adopt, 
the holographic model predicts that a normal interfer- 
ence pattern is observed at the HST (with D k 2.4m) for 
A « 1.6/im. On the other hand, if D were about 3 orders 
of magnitude larger, according to the holographic model, 
a halo structure would be observed. 



V. TESTING SPACETIME FOAM MODELS 
USING ASTRONOMICAL OBSERVATIONS 

With the above as background, we are now prepared to 
discuss the issue of testing spacetime foam models. This 
subject has been discussed previously by CNvD as well 
as [H] and [H. 

CNvD described a more flexible way to use astronom- 
ical observations, which has been summarized and ex- 
panded upon in the preceding sections. As we note in 
the Introduction, for point sources of low redshift z, it 
is not so critical which distance expressions for the point 



sources one should use. But for high redshift objects, 
it does make a difference. We have already argued that 
the appropriate cosmological distance to use is the total 
line-of-sight comoving distance (lpj given by 

D c {z) = D H f dz'/E(z'), (9) 
Jo 

where 



e(z) = ^n M (i + z)3 + n k (i + z y + n A , (10) 

with Dh = c/H being the Hubble distance, ilM,flt and 
f^A being the (fractional) density parameter associated 
with matter, curvature and the cosmological constant re- 
spectively. Consistent with the latest WMAP + CMB 
data, we will use = 0.25, fl\ — 0.75 and i\, = 0, 
and for the Hubble distance we will use Dh = 1.3 X 10 26 
meters. In terms of the comoving distance, for the vari- 
ous models of spacetime foam (parametrized by a), the 
equivalent halo size is given by 

Sip = N(l -a)l%D 1 H a I(z,a)/X (n (11) 

with 

/(*, a) = jf dz'il + z')/{E{z')) (£ dz"/E(z")^J , 

(12) 

where the factor (1 + z') in the integral corrects the 
observed wavelength A , back to the wavelength X(z') 
at redshift z'. That is A(z') = A Q /(1 + z'). We have 
used these results to produce Figure [TJ which shows 
the predictions made for three different values of a, i.e. 
a = 2/3,0.6,1/2 respectively, for the size of observed 
haloes produced by accumulated phase dispersion for a 
source at two redshifts, respectively z — 4 and z = 1. 
This graph plots wavelengths from hard X-rays (equiv- 
alent photon energy ~ 100 keV) down to radio waves. 
The labeled arrows pointing upward to the right delin- 
eate the diffraction limited response of various existing 
or proposed telescopes. That is, with the exception of 
Chandra, the arrows plot X /D, where D is the aper- 
ture/baseline length and the length of the arrow is de- 
termined by the telescope's wavelength limits. The shal- 
lower slope of Chandra 's results from its grazing-incident 
optics and its complicated point spread function. 

As Sip depends nontrivially on the redshift of the 
source, a remark on the effect of the uncertainties in the 
determination of the sources is in ordered. We have cal- 
culated the ratio S I/I for a variety of redshifts and un- 
certainties. Some sample values are given in Table [I] As 
can be seen, the effect is small, unless the uncertainty 
in the redshift, Sz, grows significantly larger than 0.1, 
as would be the case for photometric redshift determina- 
tions [m, m 

Steinbring [39| showed a plot of a somewhat similar na- 
ture based, however, on luminosity distance rather than 
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FIG. 1. The detectability of various models of foamy space- 
time with existing and planned telescopes. Diagonal tracks 
are shown for three models of foamy spacetime, namely 
a = 0.5,0.6,2/3, for z = 1 and 4 (respectively the lower 
and upper tracks for each model) and N = 1.9. Also shown 
are the observing ranges and theoretical resolution limits (i.e., 
PSF size) for a wide variety of telescopes, both current and 
planned. The tracks for each telescope represent the diffrac- 
tion limit with the exception of the dotted track for Chandra, 
whose complicated PSF is NOT diffraction limited (see sec- 
tion [V])- As such, they represent the detectability limit for 
these observations, in the case of a perfect telescope of that 
size. Importantly, note that in practical terms a given tele- 
scope may not be able to probe quite as far as pictured, as in 
the final analysis the region probed is also a function of Strehl 
ratio (see section fVT]) . 



the co-moving distance (his Figure 1). The main dif- 
ferences here are corrected curves for a = 2/3 and 1/2, 
since [3j| also did not show the curves yielded by CNvD. 

We note that if the arrow representing a telescope's 
diffraction limited response lies below the halo size curve 
for a given a, that model may be excluded by observa- 
tions. For example, as mentioned previously the random 
walk model (a = 1/2) is clearly excluded by HST obser- 
vations, as well as radio observations by the NRAO Very 
Long Baseline Array (VLBA) . The latter surveys find nu- 
merous examples of unresolved sources at cosmological 
rcdshifts [IH, [l3| . We also note that HST observations 
at short wavelengths (~ 10 -6 — 10 -7 m) appear to be 
capable of testing the holographic model (a — 2/3) and 
certainly may exclude a = 0.6. 

Before delving further into the observational record to 
discuss what models of spacetime foam are excluded by 



the current data, it is useful to make a few comments. 
First of all, what types of observations and objects would 
be useful for constraining or rejecting spacetime foam 
models? In principle, all one needs are high-resolution 
imaging observations of a source that is unresolved and 
at high redshift. However, as noted previously, for a va- 
riety of practical reasons the source must be reasonably 
bright. This is for two reasons. First, one needs to be 

- able to reliably distinguish a halo produced by disper- 
sions induced by quantum fluctuations in spacetime, from 
a faint "halo" that might be due to intrinsic structure in 
the source. For example, in the optical, every bright ex- 
tragalactic point source, be it a supernova, QSO, GRB 

1 or otherwise, is resident in a galaxy, which could very 
easily mimic or disguise a halo. Thus one wants to be 
able to achieve large dynamic ranges (at least hundreds 

- or thousands to one) on the image with a reasonable 
_ amount of integration time. This host galaxy can also 

scatter light from the AGN, so that there may still be 
some halo left in the image (e.g., [III). However, the 
" scattered light is polarized, so that there is a clear way 

- to discriminate the effects of spacetime foam from scat- 
tering effects. Secondly, a number of observing methods, 
particularly optical interferometry, require much brighter 
sources in order to be able to find fringes in real time. 
For example, the design of the VLTI requires sources in 
the IR K-band brighter than 14th magnitude. Folding 
these constraints in, for non-interferometric observations 
in the optical, any bright point source will do, be it a QSO 
or GRB. Unfortunately, SN la or SN II cannot be used, 
as they are comparable in luminosity to the host galaxy 
so they will not allow large enough signal to noise to 
be reached (i.e., we will not be able to distinguish struc- 
ture in the host galaxy from broadening due to spacetime 
foam). However, for optical interferometers such as the 
VLTI or (in the future) SIM or Darwin/TPF, only the 
brightest Blazars observed in outburst, or the prompt 
flash of a powerful GRB will do. 

In the radio, any bright QSO or BL Lac object will suf- 
fice, and fortunately there are large numbers of these. A 
large number have been found to have unresolved sources, 
as seen with both the NRAO VLA and VLBA [H El]. 
These observations confirm at radio wavelengths that the 
random walk model can be ruled out. None of these 
observations have sufficient resolution to test the holo- 
graphic model, however. 

Ideally, one wants to observe with any telescope that 
will yield the smallest X/D ratio as the diffraction limit 
would then dictate that one could set the smallest limit 
on the size of any halo corresponding to the wavelength 
at which a given spacetime foam model would predict 
the largest halo size. However, for one of the telescopes 
shown on Figure 1, namely the Chandra X-ray Observa- 
tory, this is not the case. Unlike the radio, IR and op- 
tical telescopes shown there, Chandra's optical design is 
governed by the short wavelengths (comparable to inter- 
atom spacings) of X-ray photons, which make ordinary 
normal-incidence optics impractical. Instead, Chandra 
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uses grazing incidence mirrors in a nested [H, |4{| de- 
sign. In this design, it is the nesting of the mirrors 
that focusses the light on the imaging plane. As a re- 
sult, the PSF is not diffraction limited, and the telescope 
cannot sense the wavefront coherently. We have plotted 
Chandra's angular resolution in Figure 1 using param- 
eters from the Chandra Proposer's Observatory Guide, 
using values corresponding to the 75% level in the en- 
circled energy function (see specifically Figures 4.6 and 
4.7). However, in Figure 1 we show the Chandra line as 
dotted. This is because its inability to sense the wave- 
front coherently makes it unsuitable for testing models 
of spacetime foam because the blurring of the images is 
a product of interference effects rather than photon scat- 
tering, an effect which non-diffraction-limited optics arc 
not sensitive to (see section IIV[) . 

Of course one would also like to be able to obtain con- 
straints in as many bands as possible, since finding the 
halo size as a function of frequency, in principle, allows 
one to search for the characteristic A -1 dependence in the 
size of the spacetime foam scattering disk (discussed in 
section ITV]) . Constraints obtained in only one band are 
susceptible to bias by structure intrinsic to the source 
(host galaxy or extended emission line regions, for exam- 
ple) or to the PSF of that band 

With the above provisos, then, we can now discuss the 
constraints that can currently be put on spacetime foam 
models. It is apparent that with current telescopes, all 
spacetime models with a < 0.6 are firmly excluded by 
the combined archives of telescopes in the optical, radio 
and X-ray where unresolved sources have, at some level 
(see below) been observed at every redshift up to at least 
z = 5. With HST observations, this constraint can be 
significantly tightened. Thus the random walk (a = 1/2) 
model is securely excluded at very high confidence, as 
in fact are all models with a < 0.65. The Holographic 
(a = 2/3) model is very close to being tested and will 
enter that regime once the VLTI is fully operational, and 
may be tested by HST in the ultraviolet. 



VI. TESTING SPACETIME FOAM MODELS - 
THE HST QUASAR ARCHIVE 

Because the Strehl ratio is defined as the ratio of the 
observed image peak to the peak diffraction spike of an 
unabberated telescope, it provides a useful way to assess 
the effects of spacetime foam. As discussed above, space- 
time foam causes wavefront "errors" in the images of dis- 
tant point sources. The effect is similar to atmospheric 
"seeing" and, in fact, we follow the general analysis of 
Sandler et al. 36] for the degradation of Strehl ratios 
resulting from wavefront errors. Since spacetime foam 
induced wavefront errors are not spatially correlated, a 
simple approximation for the degraded Strehl ratio is [36| 

Sobs = S M exp[— (o-j + cr|)] (13) 



where Sm < 1 represents a degradation of the observed 
Strehl ratio due to masking effects (which may or may 
not be spatially correlated as discussed above), 07 rep- 
resents uncorrelatcd wavefront errors induced by the in- 
strumentation (i.e., telescope plus instruments) and 0$ 
represents uncorrelated wavefront "errors" induced by 
spacetime foam. Both of these dispersions are expressed 
in units of the telescope's diffraction limit, X/D. For 
convenience we then define the spacetime foam degraded 
Strehl S<f, as: 

= exp(-al) (14) 

where cr^ is simply the angular dispersion of the wave- 
front vector due to external effects (i.e., Sip in Eq. (JTTJ)) 
divided by X/D, 

0-< f> =N(l-a)l%D 1 I f a I{z,a)D/X 2 o . (15) 

(S$ ~ 1, of course, corresponds to no image degradation.) 
This approximation, of course, breaks down when <tj, ~ 1, 
i.e., when the wave front angular dispersion is compara- 
ble to the telescope's resolution. A fully parametrized 
version of the resultant Strehl ratio then is 

= exp (~A^ 2 (1 - a) 2 l 2 P a Df- a) I 2 {z, a)D 2 /X^j . 

A few comments concerning parameter sensitivity of 
these models for Strehl ratio are in order: 

1. The model Strehl ratio is strongly sensitive to wave- 
length. 

2. Although we have managed (section |ll| to narrow 
the range of possible values for the constant, N, in 
the holographic model the remaining uncertainty in 
this parameter (0.97 < N < 1.9) leads to a factor of 
3.4 in the argument of the Strehl ratio exponential. 

3. The instrumental contribution, S(o~i) can be mea- 
sured using a true point source, for example a star. 

4. The model Strehl ratio is also exquisitely sensitive 
to the value of a. 

Figure [2] illustrates the sensitivity to degradation of S<p 
from spacetime foam as a function of redshift and of 
wavelength for a values of 2/3 (in black) and 0.655 (in 
grey) - assuming N = 1.9 in each case. In the left-hand 
figure we illustrate the change in expected for a point 
source with varying distance for three fixed wavelength 
values. In the right-hand figure we show the change in 
expected for a point source with varying wavelength 
at three fixed distances. We can see that HST is just 
able to probe the a = 2/3 case at the very shortest 
(UV) wavelengths in its spectral range. We note that 
N = 0.97 shifts the curves in the left-hand figure to the 
right (i.e., increasing redshift). We also note that shorter 
wavelength is a greater imperative than increased dis- 
tance for improving sensitivity to a. However, optical 
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measurements may rule out lower values of a. Clearly 
the realm for testing the holographic model lies in red- 
shifts z 1. Notice, also, that even a modestly smaller 
value of a = .64 clearly can be tested at redshifts more 
than an order of magnitude smaller. 

To test the spacetime foam models using 
Strehl ratios, we explored the HST archive 
( |http: / / archive.stsci.edu /hst / search. php) for obser- 
vations of high-redshift quasars, using StarView to 
cross reference suitable catalogues with the archive 
database. We explored the Sloan Digital Sky Survey 
(SDSS) DR5 quasar catalogue 13811 a nd a catalogue of 
the highest redshift SDSS quasars [33| in addition to the 
quasar catalogue of [44[ - hereafter VCV. In Table |TT] 
we present notable high-redshift quasars observed in the 
optical, NIR and UV using HST . We prioritize those 
observations that maximize the value of I{z,a) / \ D (see 
Eq. G2J). 

The current record for the highest redshift quasar is 
held by the CFHT quasar survey: CFHQS J2329-0301 
at z = 6.43. This has not yet been observed with HST . 
The highest redshift quasars with HST pointings all come 
from the SDSS, including SDSS J1030+0524 at z = 6.28. 
and SDSS J1148+5251 at z = 6.42 The former was ob- 
served, along with other high-z SDSS quasars, using the 
Advanced Camera for Surveys (ACS) High Resolution 
Channel (HRC) at A G ~ 850nm by (33|. Blurring of these 
images has been studied inconclusively [Hf. However, 
the observations are all at red wavelengths (F775W and 
F850LP), which impinges on their sensitivity to space- 
time foam, especially the a = 2/2, case (see Eq. HH and 
Figure [2]) ■ The latter was observed through two linear 
ramp filters by [46] at - 720 and - 900 nm, and by 
in F850LP and F775W. Among other notable examples is 
the HDF-S QSO at z = 2.24 37] , which was observed as 
part of the multi-wavelength coverage of the deep fields 
using the Space Telescope Imaging Spectrograph (STIS) 
in the ultraviolet. While offering the big advantage of 
shorter wavelengths, the STIS data is unsuitable for the 
present purpose due to the extreme distortion of the PSF, 
which appears triangular with no discernible Airy ring. 

We were unable to find any short wavelength (A < 
500 nm) images of high redshift (z > 4) quasars by cross- 
referencing with SDSS or VCV quasar catalogues. These 
short wavelength observations are critical because they 
represent the only range of the spectrum for which Fig- 
ure [1] indicates we may be able to probe the a = 2/3 
case. 

We turned next to the Hubble Ultra-Deep Field 
(UDF). The UDF Q is the deepest image taken of the 
sky to date. The 10 6 s long exposure was taken with 
the ACS using a multi-angle dither pattern to smooth 
out the PSF. A number of intermediate-to-high redshift 
quasars have since been identified in this field, with red- 
shifts provided by the GRAPES projec t [3~i| . A list of 
quasars in the UDF was published by [51(. Four such 
quasars were chosen as suitable candidates due to their 
relative isolation in the field (see Table HI)) . While we ac- 



TABLE I. Sample Uncertainties in I 



z 


8z 


a SI(z, 


a)/I(z,a) 


1.00 


0.01 


0.50 


0.0104 


1.00 


0.01 


0.60 


0.0101 


1.00 


0.01 


0.67 


0.0099 


5.00 


0.01 


0.50 


0.0018 


5.00 


0.01 


0.60 


0.0013 


5.00 


0.01 


0.67 


0.0011 


5.00 


0.05 


0.50 


0.0090 


5.00 


0.05 


0.60 


0.0067 


5.00 


0.05 


0.67 


0.0055 



knowledge the UDF data is not ideal for this task, since 
the sources will be inherently faint, and the data is as- 
sembled in a complicated stacking procedure which may 
introduce smoothing, it is only non-ideal in that it weak- 
ens our possibility of testing our null hypothesis: any 
smoothing effect serves to dilute the Strehl ratio further. 
It is clear from the table that the advantage of using 
shorter wavelength B-band data more than compensates 
for the lower redshift of the sources when compared to 
the SDSS quasars examined by [3!| . Images of each UDF 
quasar in each band (B, V, i and z) are presented in Fig- 
ure [3] It is noticeable from the original UDF images, 
that HUDF-QSO 6732 appears somewhat blurred in the 
B filter (F435W). 

As mentioned above, the Strehl ratio is defined as the 
ratio of the image peak to that of the diffraction spike 
of an unabberated telescope. The unaberrated PSF was 
calculated assuming a primary mirror radius of 1.2 m 
with an inner radius of 0.35 m, convolved with the filter- 
instrument response function at a resolution of lA to pro- 
duce an unaberrated PSF for each filter. These are shown 
along with stellar PSF's for each band in Figure [H How- 
ever, as discussed by [Hj], HST has a complicated non 
azimuthally symmetric structure to its PSF due to the 
support structure of the telescope itself. This is partially 
mitigated by the UDF observing strategy which smoothes 
the PSF and makes it rounder. Furthermore, the process- 
ing that goes into the UDF removes the distortion present 
in raw images. Note that the SNR of the quasars is too 
low to be able to obtain reliable photometry on the indi- 
vidual raw frames, which could be accurately compared 
to Tinytim models, as was done by [39]. We found it 
futile to model the UDF PSF using TinyTim due to the 
complicated stacking process ("drizzling"). Instead we 
resort to stellar PSF stars to produce "corrected" Strehl 
ratios (after Steinbring 2007). While we cannot control 
the SED of the stars, we can perform repeat experiments 
using multiple stars from across the field. 

We determined the centroid of each image and 
coaligned quasar and stellar PSF images to the centre 
of the nearest pixel. We measured the Strehl ratio for 
each quasar and stellar PSF image, with respect to the 
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FIG. 2. Expected Strehl ratio with redshift (left) and wavelength (right) for two spacetime foam cases: a — 2/3, N = 1.9 
(black lines); a — 0.655, N = 1.9 (grey lines). In the left-hand figure we illustrate the change in Strehl ratio expected for 
a point source with varying distance, z, for three fixed observed wavelengths (A = 200, 300 and 400 nm as marked). In the 
right-hand figure we show the change in Strehl ratio expected for a point source at various three fixed distances: z — 6 (solid 
line), z = 4 (dashed) and z = 2 (dot dashed). 




FIG. 3. The four HUDF quasars studied in this paper, from left to right: 9397 (z = 1.2); 6732 (z = 3.2); 9487 (z = 4.1); 4120 
(z = 2.1). Descending rows show the F435W (B), F606W (V), F775W (i) and F850LP (z) images. 



ideal, azimuthally symmetric PSF - see Table IIIII We 
then divide the quasar Strehl ratio by each Stellar PSF 
Strehl ratio to produce a corrected Strehl, Sqso/ Sstar, 
after [39| . This is found to be equivalent to measuring the 
Strehl ratio with respect to the stellar PSF, SQso-star 
- see table IIVI - which is as expected for well-centered 
images. Referring to Eq. Q21 we get 

Sstar = exp[-o-f], (17) 

so 

Sgso/Sstar = Sm exp[—a^] = SmS<p. (18) 



We find that the stars produce quite stable Strehl ra- 
tios in each filter, across the field (varying at the < 5% 
level) , while the quasars are found to vary in Strehl ratio 
considerably in a given filter. We measured all Strehl ra- 
tios using a range of photometric radii, from just 9 Airy, 
up to Of 5, and found results to be consistent within the 
errors at each photometric radius. We adopt 0.'25 as our 
photometric radius. Errors are dominated by systemat- 
ics, estimated conservatively at ~ ±10% as may be in- 
ferred from Figure [S] where a few points have Sqso / Sgtai 
ratios ~ 1.1. Repeat measurements at differing photo- 
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TABLE II. Quasar pointings with high I(z,a)/\ (see Eq. lll[) values in the HST archive. I(z, a)/A is a measure of the number 
of wavelengths travelled by the light detected at wavelength A . All calculations in this table assume a = 2/3. 

Name z Instrument Filter A G I(z,a) I(z 1 a)/X 











H 




[Mm- 1 ] 


SDSS J1044-0125(33. 39] 


5.80 


ACS/HRC 


F850LP 


0.91 


6.98 


7.71 


SDSS J-0836+0054f33. 391 


5.82 


ACS/HRC 


F850LP 


0.91 


6.99 


7.72 


SDSS J-1306+0356[33. 391 


5.99 


ACS/HRC 


F850LP 


0.91 


7.06 


7.81 


SDSS J-1030+0524[33, 391 


6.28 


ACS/HRC 


F850LP 


0.91 


7.19 


7.95 


SDSS J-0913+5919134, 391 


5.11 


ACS/HRC 


F775W 


0.77 


6.62 


8.60 


SDSS J-2228-0757f34. 391 


5.14 


ACS/HRC 


F775W 


0.77 


6.64 


8.62 


SDSS J-1208+0010f34, 391 


5.27 


ACS/HRC 


F775W 


0.77 


6.71 


8.71 


SDSS J-0231-0728f34. 391 


5.41 


ACS/HRC 


F775W 


0.77 


6.78 


8.81 


SDSS J1148+5251[46j 
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metric radii gave far lower variation ^ 1%. 

The Strehl ratios of each quasar from the UDF with 
respect to the UDF stellar PSF's are plotted in Figure [SJ 
against the theoretical Strehl ratio obtained for various 
values of a and N. Any measured degradation above the 
1:1 line (shaded region) rules out that spacetime foam 
model, since the quasars are not sufficiently blurry. That 
is, the image is degraded less than predicted by the the- 
oretical model. Points below the 1:1 line may indicate 
structure in the sources, i.e., the blurriness is above what 
we would expect for a given spacetime foam model and, 
thus, does not test the theory. 

Some blurring is seen, consistent with the level ex- 
pected if a = 0.655. All models with a < 0.65 are ex- 
cluded by the UDF observations. 

It is straightforward to disprove a given model, by find- 
ing an unblurred quasar just a little more distant or at 



shorter wavelengths (higher in the plot). Points lower 
in the plot may imply structure in the source. However, 
with the data we currently have, the alternate explana- 
tion, namely the combination of limited Strehl ratio and 
the signal to noise of the UDF quasar observations, is 
equally plausible. We are clearly not yet probing the 
a = 2/3 scenario. Larger statistical samples eventually 
may allow us to disentangle blurring due to structure 
from that due to the effects of spacetime foam. 



VII. DISCUSSION AND CONCLUSION 

Using the UDF data, we see evidence for an increased 
blurring (decreased optically corrected Strehl ratio) with 
increasing I(z,a)/X . This is likely due to masking or 
source structure in the majority of cases, but we note 
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FIG. 4. Top, from left-to-right: HUDF-QSO 6732 in F435W (B), F606W (V), F775W (i) and F850LP (z) bands. Middle: The 
ideal model PSF (1.2 m radius mirror with central hole of 0.35 m radius) for same filters. Bottom: Stellar PSF's for the same 
filters. 



TABLE III. Strehl ratios for UDF quasars (with identifier 
from [5j|) and stars (with UDF coordinate positions) in each 
filter, with respect to the idealized, azimuthally symmetric 
PSF. Unmeasurable values occur where the image is satu- 
rated. 



Source 


F435W 


F606W 


F775W 


F850LP 


QSO 9397 


0.31 


0.40 


0.60 


0.60 


QSO 6732 


0.24 


0.37 


0.56 


0.52 


QSO 9487 




0.17 


0.28 


0.29 


QSO 4120 


0.25 


0.35 


0.44 


0.56 


Star A (3026,3408) 


0.27 






0.64 


Star B (5692,4830) 


0.28 






0.63 


Star C (5006,2496) 


0.28 


0.40 


0.59 


0.62 


Star D (8459,6278) 


0.28 


0.34 




0.62 


Star E (1231,5437) 


0.29 


0.36 


0.58 


0.63 



TABLE IV. Strehl ratios for the UDF quasars (identifiers 
from [5l[) with respect to stellar PSF's, Sajso-Star, using a 
photometric aperture of 0.'25. 



Source 


F435W 


F606W 


F775W 


F850LP 


9397 


1.10 ± 0.1 


1.02 ± 0.1 


1.02 ± 0.1 


0.98 ± 0.1 


6732 


0.88 ± 0.1 


0.94 ± 0.1 


0.95 ± 0.1 


0.86 ± 0.1 


9487 


0.00 ± 0.5 


0.42 ± 0.1 


0.48 ± 0.1 


0.47 ± 0.1 


4120 


0.90 ± 0.1 


0.88 ± 0.1 


0.75 ± 0.1 


0.92 ± 0.1 



that the blurring is at a level consistent with a spacetime 
foam model with a = 0.655. All models with a < 0.65 
are excluded by the UDF observations. 



We are not yet probing the a — 2/3 regime, however; 
the UDF data is non-optimal for this purpose as they are 
not at a wavelength short enough to exclude the holo- 
graphic model. Indeed, we see obvious blurring in the 
quasar images, in some cases due either to structure or 
masking. Nevertheless, the UDF observations provide 
the strongest constraints that currently exist on space- 
time foam. 

We have shown that the current HST archive does 
not contain observations that can more tightly constrain 
spacetime foam models. To do this, we would need ob- 
servations in the ultraviolet. While one such observation 
currently exists, it was made by the STIS instrument, 
which has a triangular PSF with no Airy ring, making 
the analysis we detail in section IVT1 impossible . The mea- 
surement can be done cither by the new WFC3 instru- 
ment, which can access the very edge of the a = 2/3- 
sensitive region of the parameter space at the blue end of 
its spectral range (using the F225W, F275W or F300X fil- 
ters), or alternately by the ACS SBC, which can probe to 
100 nm. While both instruments under-sample the PSF 
in the ultraviolet, each has its advantages. The WFC3 
is considerably more sensitive, so that higher dynamic 
ranges can be probed in a shorter exposure time. How- 
ever, at the short wavelengths probed by the SBC the 
anticipated blurring is so strong that the under-sampling 
may not matter. 

The two main difficulties in probing high redshift 
quasars at such short wavelengths are scattering by their 
extended haloes, and absorption by intervening HI and 
Hell. The former is the most likely explanation for any 
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FIG. 5. Measured Strehl ratio S for the four UDF quasars (9397: circles & solid line, 6732: up-pointing triangles & dotted line, 
9487: squares & dot-dashed line, 4120: down-pointing triangles & dashed line) against theoretical Strehl ratio (Eq. I16p . for 
various spacetime foam models. The shaded area indicates S > S<p: Such observations exclude a given spacetime foam model, 
(a, N). Crosses indicate data points from [3S|. We exclude all models with a < 0.65 since we clearly see quasars that are less 
blurred than is predicted by such models (top). We do in fact see a small amount of blurring in our UDF quasars, consistent 
with the magnitude of blurring expected for a = 0.655 (middle). No degradation should be seen for a — 2/3 (bottom) for the 
wavelengths and distances probed so far by the HST archive (A > 400 nm). 
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observed structure in the quasars, should it be observed, 
but can be distinguished from the effects of spacetime 
foam using polarimetric measurements (e.g., [U 
Furthermore, we do not expect scattering to dominate 
the flux for type 1 quasars, even at wavelengths ~ 1215A, 
since [43[ found a wavelength independent scattering law 
for obscured quasars. The latter is the more signifi- 
cant problem, as we are already photon starved at rest- 
frame wavelengths shortward of Ly alpha, and at red- 
shifts > 4.3 both HI and Hell absorption will affect the 
opacity at wavelengths < 1215A in the quasar frame 
[171 Il9| . Nevertheless, in the absence of diffraction- 
limited X-ray imaging in the foreseeable future, rest- 
frame far-UV imaging represents the best hope for prob- 
ing the a — 2/3 regime, although we admit we would 
be lucky to find a quasar sightline that unambiguously 
disproved the a = 2/3 case by detecting no blurring or 
structure. In the rather more likely event of structure be- 
ing detected, it would take a series of deep and expensive 
observations to establish that the blurring or structure 
was not due to scattering. The best sources to inves- 
tigate are the high resdshift quasars from the Richards 
(SDSS) catalogue that are known to be isolated and to 
show no structure in their z or I band images. All the 
UDF quasars show companions or similar nearby struc- 
ture. 

As already noted at the end of section [V] one of the 
most exciting prospects is that the fully operational VLTI 
will allow us to test the Holographic model. In fact, as 
shown by Figure [TJ due to the combination of its base- 
line size plus wavelength, the VLTI emerges as the best 
overall current or planned instrument for the testing of 
spacetime foam models. The instrumental setup to use 
would be cither MIDI (at wavelengths longer than 5 /jm) 
or PRIMA (at shorter wavelengths). This would be a 
rather difficult observation, because currently, VLT in- 
terferometry is limited to the very brightest sources: in 
K band, the practical limit is K — 10, although that can 
be extended down to K = 15 if a bright star of K = 10 
or brighter is very close to (within a few arcseconds) the 
quasar (F. Dclplancke, priv. comm.). The presence of 
such a bright point source allows one to 'phase-reference' 
the observation so that the visibility of the fringes can 
be tracked with high precision. This is a significant diffi- 
culty, as the number of quasars brighter than K = 10 is 
very small - in fact, a search in SIMBAD revealed only 
one such object, namely 3C 273, which is at the modest 
redshift of z = 0.158. If one looks for bright stars within 
40 arcsec, one finds a few other quasars with bright stars 
nearby, specifically PKS 0435-300, HS 1227+4641, OK 
568, AH 26 and 2GZ J011339-3343. All of these are at 
the lower end of the useful range as far as the quasar 



brightness is concerned, and also, two or three are too 
far north for VLTI observations (HS 1227+4641, OK 568, 
possibly AH 26). 

The initial conclusion would then be that this is a dif- 
ficult observation to do and not many sources are avail- 
able to do it. However, one can use targets of opportu- 
nity to explore further. For example, bright blazars in 
outburst can reach up to K ~ 12 and in extreme cases 
K ~ 11, making them possibly useful for these obser- 
vations, particularly if the future holds an increase in 
sensitivity. Another possibility is to use the afterglows 
of bright gamma-ray bursts, which have now in at least 
one case gotten considerably brighter than 10th magni- 
tude in K. Examples here include GRB 990123, which 
reached V ~ 9 in its prompt phase (Galama et al. 1999), 
and GRB 060607A (Ziaeepour et al. 2008), which was 
B ~ 13 in its prompt phase. However, considering that 
GRBs decline precipitously in brightness within minutes, 
their utility for VLTI observations is likely quite limited. 

We conclude with a succinct summary. In this pa- 
per we elaborate on an earlier proposal (CNvD, 2006) to 
detect spacetime foam by looking for seeing disks in the 
images of distant quasars and AGNs. Assuming isotropic 
fluctuations, we argue that spacetime in effect creates a 
seeing disk whose angular diameter is ~ 6ijj and that the 
effect is expected to show up in interferometric fringe 
pattern as decreasing fringe visibility and reduction in 
Strehl ratios as soon as Sip ~ X/D. A model of space- 
time foam is disproved if images of a distant source do 
not exhibit the blurring predicted by the model. Thus 
far, images of high- redshift quasars from HUDF provide 
the most stringent tests of spaceteime foam models. The 
random-walk model is convincingly ruled out while the 
holographic model is still viable. We anticipate that the 
latter model can be tested when the VLTI is fully oper- 
ational; it is also possible that observations with HST in 
the ultraviolet may accomplish this goal. 
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